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(g) Verfahreh zur Mpmentarifrequenz-Detektion 

(§) Die Erf indung betrifft ein digital arbeitendes Verfahren zur 
Momentanf requenz-Detektion insbesondere fur hohe Fre- 
quenzen(> 100 MHz) und kiirze Me&dauer (< 100 ns), die 
insbesondere bei FSK- und/oder PSK-Verfahren aber auch 
ailgemein bei Kommunikations- und/oder Radarverf ahren 
auftreten. Dabei werden analoge Etngangssignale zunachst 
phasenquantisiert und anschlie&end mit Hirfe einer Mittel- 
wertsbildung in zugehdrige dtgttale Frequenzworte umge- 
wandelt. pamit ist insbesondere eine sehr hohe Frequenz- 
auflosung mit einem sehr geringen und kostengunstigen 
Schaltungsaufwand erreichbar. 
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. Beschreibung t ( " ' ■ ' > y 

Die Erfindung betrifft ein Verfahren zur Momentanfrequenz-Detektion nach dem Oberbegriff des Patentan- 
spruchs 1. Insbesondere in der Radar- und/oder Kommunikationstechnik, vor allem bei der Demodulation von 
5 sogenannten Frequency Shift Keying (FSK> und/oder sogenannten Phase Shift Keying (PSK^Signalen, ist es 
erforderlich,die Mpmentanfrequenz eines analogeh elektrischen Signals zu bestimmen. 

Ein Verfahren besteht darin. mit Hilfe einer sogenannten Torschaltung (Start-Stop-Schaltung) die gesuchte 
Momentanfrequenz zu bestimmen. Dabei wird bei dem analogen (Wechselspannungs-) Signal z. B. zwischen 
zwei Nulldurchgangen, z. B. jeweils zwischen den ansteigenden Flanken, die Zeit gemessen und daraus die 
io Momentanfrequenz bestimmt. , 

Diese Zeiterfassung erfolgt im allgemeinen in digitaler Form mit Hilfe eines Taktsignales, dessen bekannte 
Pulswiederholfrequenz wesentlich hoher ist als die hochste Periodendauer des zu rhessenden Signals. Es wird 
dann lediglich die Anzahl der Takte zwischen den erwahnten Nulldurchgangen des Signals gezahlt. Aus dieser 
Anzahl und der bekannten Pulswiederholfrequenz kann die gesuchte Frequenz bestimmt werden. Diese Anzahl 
15 kann aber auch als ein digitales Frequenzwort betrachtet werden, welches der gesuchten Frequenz zugeordnet 
ist. Ein solches Frequenzwort kann z. B. in einer nachgeschalteten Datenverarbeitungsanlage in digitaler Form 
ausgewertet werden. 

Ein erster Nachteil eines solchen Verfahrens besteht darin, daB die Pulswiederholfrequenz aufgrund des 
Abtasttheorerris immer wesentlich groBer als die zu messende Frequenz sein muB. Die hdchstmogliche Pulswie- 
20 derholfrequenz ist aber begrenzt, insbesondere durch die Elektronenlaufzeiteri in den verwendeten elektrbni- 
schen Bauelementen und/oder Baugruppen. 

Ein zweiter Nachteil besteht darin, daB Quantisierungsfehler entstehen, die insbesondere bei hohen (Signal- 
)Frequenzen storend sind. Diese puantisierungsfehier beruhen z. B. darauf, daB das zu bestimmende Signal und 
das Taktsignal im allgemeinen unkorreliert sind Zu diesem systematischen Quantisierungsfehler kommen dann 
25 insbesondere bei hohen Frequenzen noch weitere unsystematische Fehler, die z. B;, auf Rausch- unoVoder^ 
Jittereffekten realer Schaltungsanordnungen beruhen. ; 

Der Erfindung liegt die Aufgabe zugrunde, ein gattungsgemaBes Verfahren anzugeben, mit dem in zuverlassi- 
^ger und kostengunstiger Weise auch bei hohen Frequenzen sowie wahrend einer moglichst kurzen MeBdauer 
eine Frequenzbestimmung moglich wird. 
30 Diese Aufgabe wird gelost durch die im kennzeichnenden Teil des Patentanspmchs 1 angegebenen Merkmale. 
Vorteilhafte Ausgestaltungen und/oder Wei terbildungen sind den Unteranspruchen entnehmbar. 

Ein erster Vorteil der Erfindung besteht darin* daB bei dem zu messenden Signal zunachst eine Phasenquanti- 
sierung vorgenommen wird und die dabei erzeugten digitalen Phasenworte anschlieBend mit Hilfe einer digita- 
len Mi ttel wertbildung ausgewertet werden: Dadurch wird eine hohe MeBgenauigkeit erreicht. . . 
35 Ein zweiter Vorteil besteht darin, daB eine Signalverarbeitung in einer parallelen Grund-Schaltungsstruktur 
moglich wird. Dadurch ist eine extrem schnelle Signalverarbeitung moglich, z. B. konnen mit derzeit erhaltlichen 
Halbleiterschaltungen analoge Signale mit einer Frequenz groBer 100 MHz und einer MeBdauer kleiner 100 ns 
zuverlassig detektiert werden. 

Ein dritter Vorteil besteht darin, daB fur die Mittelwertbildung keihe Koeffizienten-Multiplizierer benotigt 
40 werden, 

Ein vierter Vorteil besteht darin, daB eine Verringerung der Abtastwerte moglich wird ohne daB eine 
Verringerung der Frequenzaufldsung und/oder MeBgenauigkeit erfolgt. Dadurch sind auch hochste iFrequenzen 
mit relativ langsamen Halbleiterbauelementen meBbar. 

Ein funfter Vorteil besteht darin, daB eine Mittelwertbildung mit einer Kaskadierung von Teilfiltern moglich 
45 wird, wodurch eine exponentielle Steigerung der stochastischen MeBgenauigkeit bei gleichzeitiger Reduzierung 
der Speicherzellen — und/oder Add ierstufenanzahl erreicht wird. 

Weitere Vorteile sind der nachfolgenden Beschreibiing entnehmbar. 

Die Erfindung wird im folgenden anhand von Ausfuhrungsbeispielen und unter Bezug auf schematisch 
dargestellte Zeichnungen naher erlautert. Es zeigen 
50 Fig. 1 bis Fig, 7 schematisch dargestellte Diagramme zur Erlauterung der Erfindung. 

Die Erfindung beruht auf der direkten digitalen Detektion (DDD) von Momentanfrequenzen, was nachfol- 
gend naher erlautert wird. 

Das Ziel einer digitalen Momentanfrequenz-Detektion ist es, bei moglichst kurzer MeBdauer x (typisch < 1 00 
ns) eine unbekannte Frequenz f x eines von der Zeit t abhangigen sinusformigen Signals 

55 . • • 

s(t) == A x - sin(27tf x * t -f (p x ) mit 0 < <p x ^ 2n 
bzw. 

60 s(t) .« A x - sin((p(t))mit0<<p(t)<2n, (1) 

mit entsprechender Genauigkeit bzw. Auflosung q zu bestimmen. D.h. die analoge Signalfrequenz f x wird durch 
diskrete Frequenzwerte F x aus einem endlichen Wertevorrat Qf dargestellt Dabei gehort zu der Frequenz f x des 
sinusformigen Signals s(t) eine Periodendauer T x gemaB der Formel 
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Alternattv dazu kann die Momentanfrequenz'f x = f(t) von dem in GL 1 (Gleichung 1) genannten Phasengang 
<p(t)des MeBsignales s(t) abgeleitet werden. 

Mathematisch steiit die Momentanfrequenz f(t) eines Signals die zeitliche Ableitung seines Phasenwinkelyer- 
laufs cp(t), also^die Winkelgeschwindigkeit to, dar. Dam it kann als Detektionsprinzip die Bildung des Differential- )0 
quotienten im TiefpaBbereich 

f_(t)= —_ . dyT(t) mit 9 T (t) - <P 0 + 2 * J 'tW*' (3) ■' I5 

2* dt "* . 

oder im BandpaBbereich ' ' ■ 

X d<p (t)-' to . 20 

t (t)=f„+ .— mit V (t)=2tff 0 +V 0 + 27T J f T (t)dt (4) 

27T . dt 

aus dem Phasengang qn(t) des TiefpaBsignals si<t) aufgefaBt werden. Jedoch ist fur die Frequenzrhessung im 25 
BandpaBbereich im MeBergebriis fur. f(t) noch eirie; Frequenz-Offset-Konstante f 0 hinzuzuaddieren, die die 
-Mittenfrequenz des MeB-Bandpasses Bm darstellt 

Nun kann diese MeBvorschrift, GL 3 und Gl. 4, bespnders vorteilhaft bei der direkten digitalen Detektion 
(DDD) angewandt werderu Dieses Verfahren erfordert zunachst eine Phasenquantisierung des Eingangssignales 
s(t) durch einen Phasen/Digital-Urnsetzer (PDU), der z. B. aus der deutschen Offenlegungsschrift DE- 30 
OS 38 33 005 bekannt ist Dabei werden zu dem Eingangssignal s(t) gehorende digitale Phasenworte <D(k) 
erzeugt, wobei k eine ganze Zahl ist und nachfolgend noch naher erlautert wird und wobei eine Kodierung nach 
dem Johnson-Code erfolgtl 

Fig. 1 zeigt beispielsweise ein von der Zeit t abharigiges Eingangssignal s(tX das zunachst mit Hilfe von 
beispielsweise acht ( 1 < i < 8) parallel geschalteten Phasenstellgliedern, die jeweils eine Phasenverschiebung i • 0 .35 
(entspricht i-x, wobei x die 0 entsprechende Zeitverschiebung ist) bewirken, in sogenannte 8-bh-Johnson-Se- 
quenzen sq,( t), 1 ^ i ^ 8, umgewandelt wird Erfolgt nun eine digitale Abtastung dieser 8-bit-Johnson-Sequenzen 
sqi(t), z. B r jeweils in der Mine der durch die gestrichelten Linien dargestellten Zeitintervalle, mit einem Takt T, 
dessen Pulse gezahlt werden (1 k k^ »), so entstehen nach der Abtastung Phasenworte Oi(k) im 8-bit- Johnson- 
Code. In Fig.l gehort beispielsweise zu der Phasenlage bei den Takten k, k+1, k+2 die Phasenworte 40 
O(k)=inillllV<K k +0 =01 ^ im Johnson-Code vorliegenden Phasenworte 

konnen dann in eihen fur die weitere Signalverarbeitung geeigneten Code, z. B. einen Bihar-Code umgewandelt 
werden- Besonders vorteilhaft ist eine Umwandlung in den an sich bekannten 2er-Komplement-Code t was noch 
erlautert wird. 

Fig. 2a zeigt in einem schematischen Blockbild eine Prinzipschaltung mit einem vorstehend beschriebenem 45 
Phasen/Digital-Urnsetzer (PDU) fur ein analoges Eingangssignal s(tX das als BandpaBsignal vorliegt und das 
zunachst mit einem an sich bekannten Hilbert-Transformator in Quadratursignale x(t), y(t) umgewandelt wird. 
Nach der digitalen Abtastung werden die Phasenworte 0(k) mit Hilfe eines Phaseh-FrequenzAvortwandlers, der 
noch beschrieben wird, in Frequenzworte F-r(kX welche den Frequenzen des BandpaBsignales entsprechen, 
umgewandelt. Zu diesen Frequenzworten Fr(k) wird dann noch in digitaler Form ein Frequenzwort F 0 (Fre- 50 
, quenz-Offset) addiert, welches die Frequenzlage (z. B. die Mittenfrequenz) des BandpaBsignales enthalt Nach 
dieser Addition entstehen dann Frequenzworte F(kX welche absolute Frequenzangaben beinhalten. Liegt das 
Eingangssignal s(t) als Tief paBsignal vor, so wird die Addition des Frequenzwortes F 0 nicht bendtigt 

Zur Momentanf requenz-Bestimmung nach dem DDD- Verfahren ist die digitale Ausfuhrung des Dif ferential- 
quotienten des Phasenverlaufs <p(t) nach GL 3 oder 4 notig. Dies geschieht anhand des vom PDU gelieferten, 55 
digitalen Phasenwinkels <D(k) (Phasenworte). Damit kann der Differentialquotient fur f(t) durch den digitalen 
Differenzquotienten 

1 0(Jc)-0(k-n) 1 A0 k 60 

f (k)= . = - — . - — (5) 

27T T(k)-T(k-n) 27T AT 

angenahert werden. Diese direkte digitale Detektion nach Fig. 2 hat in Abhangigkeit der Phasenwinkelauflosung 65 
0 des PDU's, der Abtastperiodendauer T bzw. der zeitlichen Differenz AT der Phasenworte <D(k) sowie der 
verwendeten Rechenwerksparameter einen zugehorigen linearen Frequenzfehler Af bzw. eine endlich genaue 
Freqiienzauflosung q (Rechenwerksaiiflosung). 
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Deshalb ist fur die technische Anwendung das DDD-Grundprinzip nach GL 5 wegen seiner zu geringen 
Rechenwerksauflosung 

8 _ . 1 

— = 2 N -f A mit f A = — (6) 
T T 

meist so nicht anwendbar. Daher wird durch eine nachfolgende, digitate Integration bzw. Mittelung die Fre- 
io quenzauflosung q und weitere Parameter, wie z. B. der Signalstorabstand Df der Detektion, erheblich verbessert. 
Hierfiir wird 'nach Fig. 2b ein spezieller Mitteiwertbildner mit einer Abtastraten-Reduktion vom Faktor R und 
im Falle einer Bearbeitung von BandpaBsignalen ein Frequenzoffset-Addierer dem eigentlichen DDD-Detektor 
nachgeschaltet. Diese vorteilhaften Detektionsfilter fur den Mitteiwertbildner werden noch naher beschrieben 
H ier soil nun zunachst auf eine technische Besonderheit bei derdigitalen Differenzquotienten-Bildung ausden 
15 Phasenwinkeln O(k) nach GI. 5 eingegangen werden. Es ist der "Modulo-Charakter" dieser Differenzquotienten 
zu beachten. Die Division, GL 5, ist also eine Zeitdivision eines digitalen Differenzwertes. Hierbei stellt der 
Divisor die aquivalente Abtastperiodendauer AT des vorliegenden Abtastsystems dar. Daraus folgu daB nur 
digitale Differenzen aus entsprechend verzogerten Phasenworten <l>(k) in einem digitalen Rechenwerk zu bilden 
sind. D.h. an den Subtrahiererausgangen des digitalen Rechenwerks, Fig. 2b, durfen keine Phasendifferenzen > 
20 1 180° | auftreten. Damit ist also eine Modulo-Subtraktion bezQglich 360° (27t) auszufuhren. AuBerdem bedeutet 
dies. daB das Abtasttheorem ~ mindestens 2 Abtastwerte pro Phasenwinkelpenbde — dadurch erfullt ist. 
Folglich lautet jetzt die GL 5 vollstindig: 

[ t— — ] mod 271 (7) 

; 27T AT ' ' \ ■ 

Die Abkurzung "mod" stellt hierbei den Moduloperator dar, wobei allgemein fflr die Modulooperation 
30 bezuglich 2n gilt: 


<p mod 2K = <p - 2(p [int — — J (8) 

2JT 


und sofnit den Restwinkei innerhalb des 2n-Bereichs bei voriiegendem, ganzzahligen Vielfachen von 2n ergibt 

Nun ist der erwahnte 2er-Komplement-Code wegen seiner zyklischen Eigenschaften besonders vorteilhaft 
zur Darstellung der mit 27i-periodischen Phase der Modulo-Subtraktion geeignet. Unter Anwendung des 2er- 

40 Komplements in der Form nach Fig. 2c lassen sieh sowohl Freqiienzen als auch Phasendifferenzen bei Vorzei- 
chenwechsel aufeinanderfolgender Phasenwerte (z, B. Obergang von Quadrant II nach Quadrant III) einfach 
ermitteln. Nach den Rechenregeln des 2er-Komp!ements wird die Differenz A—B dadurch gebildet, daB man 
samtliche Bits vom Wort B invertiert und das so gebildete Wort zum Wort A addiert sowie der Summe 
auBerdem noch den Wert + 1 (Carry) hinzufugt. Damit sind beim DDD-Verfahren keinerlei Korrekturen notig 

45 urn irgendwelche unerwunschten Phasenfehler bei Abarbeiten des Rechenwerksalgorithmusses, GL 7, zu kom- 
pensieren. Dies ist unter anderem auch ein besonderes MerkmaJ des DDD-Verfahrens. 

Die direkte digitale Detektion (DDD) von Momentanfrequenzen, Fig. 2b, erfordert in den meisten Fallen noch 
eine digitale Nachverarbeitung des MeBwertes Fi{k) mit speziellen Mittelungsfiltern. Der Diffefenzquotient, GL 
5, des DDD-Verfahrens liefert vergleichsweise nur eine geringe Frequenzauflosung q, GL 6, die bei weitem nicht 

so ausreicht, die in der technischen Realitat auftretenden Fluktuationen (Quantisierungsfehler, Rauschen, Jitter und 
Modulationseffekte) der einzelnen Phasenworte 0(k) zu kompensieren. 

Wird zur technischen Realisierung das Prinzip von ungewichteten Mittelwertbildungen, cLh. von linearen 
Mittelungen mit den Wichtungsfaktoren +1,0 und - 1 angewandt, so ergeben sich besonders einfache Rechen- 
werksstrukturen. Natuf lich mussen diese Strukturen und dazugehorigen Rechenwerksparameter, wie z. B. die 

55 Blocklange M (Anzahl der zur Mittelung herangezogenen Phasenworte <D(k)) in Abhangigkeit der Abtastfre- 
quenz fA bzw. Abtastperiode T, der Phasenauflosung 0 fur die erforderlichen Frequenzauflosungen q Q und qf 
sowie fur die Detektionsfehler Af und 07 entsprechend optimiert werden. Auch die Rechenzeit (MeBzeit) Tm und 
Laufzeit tm der Filter auf Frequenzanderung (Modulationen) des MeBsignals s(t) einschlieBlich ihrer 3-dB-Filter- 
bandbreiten f g ist in Rechnung zu stellen. 

60 Zur Filter-Dimensionierung werden unter anderem noch folgende Parameter herangezogen: 

a) die Rechenwerks-Frequenzauflosung q 0 — kurz nur als Rechenwerksauflosung bezeichnet — f die den 
kleinsten noch darstellbaren Frequenzwert infolge des endlichen Ganzzahlwertebereichs des Rechenwerks 
angibt; 

65 b) den Effektivwert at (Standardabweichung) des absoluten FrequenzifneBfehlers Af. Daraus kann dann in 

der Regel die effektive Frequenzauflosung qr des Detektionsfilters berechnet werden; 

c) der Signalstorabstand Df, der sich am Ausgang des Filters bei idealem, rauschfreien Eingangsphasensighal 
einstellt: 
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e)die effektive Bitzahl Nr.Sie ist eine fiktive Bitzahl einer angenommenen Phasenquantisierung, die bet der 
Einfachmittelung nach Fig. 3a und 3b und gleicher Filterlange M.die gleiche effektive Frequenzaufldsung qr 
ergibt, die dem zugrundeliegenden Filter mit der tatsachlichen Qiiahtisierungsbitzahj N entspricht. 

Somit gilt fur diese Detektionsfilter unter Verwendung.der effektiven Bitzahl Nf allgemein filr ihre effektive 
Frequenzauflosurig: 

1 * - ! 

q f = - . f A . 2" N f (9) 

M 

Dagegen kann nur im Fall von nicht kaskadierten Filtern und damit ohne Abtastratenreduktion die 3-dB-Fil- 
terbandbreite mit 


M 

f„ = — • f A fur M > 2 bzw. (10) 

• M(2M+1) -2 

f- ~ : • f A fur M » 2 . . 

y (2M+1) 

angegeben werden. Zwischen der Filterbandbreite f g und der effektiven Frequenzaufldsung qr besteht ein nicht , 
ganz linearer Zusammenhanjg, der noch spater aufgezeigt wird. Mit Hilfe dieser Parameter ist es nun Ieicht 
moglich, die Detektionsfilter grundsatzlich zu klassifizieren. Dies wird in den folgenden Abschnitteri kurz 
ausgefuhrt. 

Die drei Filterversionen nach Fig. 3 stellen die einfachsten, technisch moglichen AusfQhrurigsformen dar. 
Hierbei sind die beiden ersten Versionen, MSA- und ZSA-Filter, Fig. 3a und 3b, in der Anwehdung vollstandig 
gleichwertig. Beide Realisierungen haben dieselbe Obertragungsfunktion, die eine geht aus der anderen durch 
Transponieren hervpr. Derri Modulpsubtrahierer, der den eigentlichen Differenzquotienten AO(k)yGL 7, bildet, , 
ist zur Mittelung der Frequenzworte 

1 . M-l • 

F(Jc) = • 2 A0(k-i) mit (11) 

27TM ■ ■ i=0 


A0(k) = [0(k) - 0(k-l) J mod 2tt 

jeweils ein Transversalfilter mit den insgesamt gleichen Koeffizienten von dem Wert « . 1 (deshalb in Fig. 3 nicht 
eingezeichnet) nachgeschaltet Hierbei ergibt sich eine Rechenwerksauflosung 

q Q = q f = — • f A • 2" N , .. (12) 

M 

die mit der effektiven Frequenzaufldsung qr identisch ist Gegenuber der direkten Methode, GL 6, wird der Wert 
fur die Frequenzaufldsung q 0 und qf um den Faktor M verbessert, bzw. der effektive Frequenzfehler (Standard- 
abweichung) 

a = . • mit 0 » 2JT-2 14 (13) 

um den Wert 1/ i/M verringert. Damit kann fur den Signalstdrabstand dieses Frgquenzdetektors.die Formel 
D f » 6,02(N-i) + lOlogMindB (14) 

angegeben werden. Fur die Filterlaufzeit tm, die die Halfte der MeBzeit Tm betragt, ist dann zu setzen 
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r M = = • T (15) 

2 2 

5 ■ 

/ Jedoch fur die Filterbandbreite fg gilt weiterhin die Gl. 10. 
Eine weitere Variante dieses Einzelwert-lntegrationsverfahrens ist nun nach Fig. 3c moglich, bei der M — 1 
Differenzquotienten AOi(k) aus entsprechend verzogerten 0>(k)-Worten gebildet werden, die dann anschlieBend 
nur aufzusummieren sind. 
10 Die Freqijjenzworte des DAS-Filters 

1 M-2 

F(3c) = 2 [0(k-i)-0(k-i-l)] mod 2*. (16) 

15 27T(M-1) i=0 

sind formal mit dem Ausdruck 

20 l ■ . . 

F*(k) = — [**(*) - 0*()c-M-i)] mod 27T (17) 

27T(M-1) 

25 identisch, jedoch geht hierbei der Mittelungseffekt der Mehrsummanden-Summation fur die stochastischen 
Phasenwinkelabweichungen (Phasenrauschen usw.) verloren. Deshalb darf zur Bestimmung der statistischen 
Parameter nur die GL 16 herangezogeri werdenl Damit besitzt hier ciie Rechenwerksaufldsung 

30 q_ = 2£_ = . f . 2~ N (18) 

° 2 (M-l) 

eirien halb so groBen Wert wie die effektive Frequenzauflosung qf. Dies zeigt sich auch im Effektivwert des 
35 Frequenzfehlers 

i e i 

0 ± = — — • — — (19> 

27T Je 1 TJFT 


40 


und des Signalstorabstandes 

45 Df « 3,01(2N-3) '+ 101og(M-l)indB (20) 

Eine erhebliche Verbesserung der stochastischen Parameter, so vor allem beim effektiven Frequenzfehler of, 
bei der effektiven Frequenzauflosung qf und dem Signalstdrabstand Df, erbringen die im nachfolgenden Ab- 
schnitt beschriebenen Verf ahren, sogenannte Multiblock-Integrationsverfahren. 

so Grundlage des Multiblock-Integrationsverfahrens ist das sogenannte Doppelblock-Integrationsverfahren 
oder auch "Bicentroid- Verf ahren" Gleichzeitig bei geringem technischen Aufwand erbringt es nach Fig. 4 eine 
besonders gute U nterdruckung von Rauscheffekten. Es erreicht fur nicht zu groBe Phasenwinkelwortfolgen 
der Lange M (M < 200) fast den gleichen sogenannten "Processing Gain" des bekann ten Verf ahrens nach dem 
Prinzip des kleinsten Fehlerquadrates nach GauB. Jedoch ist hierbei der Rechenaufwand im Verhaitnis zum 

55 Verfahren nach GauB vorteilhafterweise wesentlich kleiner. 

Ebenfalls sind hier wie bei Einzelwert-Integrationsverfahren, GL 11, zwei technische Ausfuhrungsformen 
moglich, die identische Obertragungseigenschaften besitzen. Das Filter mit der kanonischen Struktur, Fig. 4a, 
wird deshalb als BS A- Filter und das nicht kanonische Doppelblockfiltef, Fig. 4b, als BZA-Filter bezeichnet. Bei 
beiden Detektoren werden blockweise getrennt zwei gleich lange Wortfolgen der Lange M/2, Block A urid 

60 Block B, von Phasehwinkeln aufsummiert und dann anschlieBend die beiden Blocksummen A und B (Phasenwor- 
te OA(k) und Ob(1c)) vorietnander subtrahiert. 

Die Modulosubtraktion erfolgt hier entsprechend dem Wert M • ju Damit besitzt dieses Verfahren einen 
quadra tischen Glattungseffekt zur Gewinnung der Frequenzwerte 
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4 M-l 

F(3c) = 2 [0(k-i) - 0(k-M-i)] mod M-tt (21) 

271M 2 i=0 
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bei einer Rechenwerksaufldsung qo,die mit GL 12 der einfachen Einzelwert-Mittelung identisch ist Die besonde- 
re Wirksamkeit des Verfahrens dokumentiert sich aber in de , - ; ; 

q f = * f A • 2 " N (22) 

M 

sowie bei Effektivwert des Frequenzfehlers ■ 
10 4 

a f = — — • ■■ ■ • (23) 

27T TTF 1 M/5" 1 ; 

und dementsprechend auch bei seinem Signalstorabstand 
D f ~ 6,02(N - 3) + 30 log M in dB (24) 


Rl '= 1/ R i+l = M i (25) 
N s 

R gesamt = n R i nit 1=1,2,3, ...,^-1. 
1=1 
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Mit den folgenden Beispielen, bei denen eine Kaskadierung der bisher beschriebenen Filter fur die Mittelwert- 
bildung erfolgt, ergeben sich weitere vorteilhafte Verbesserungen der Parameter bei gieichzeitiger drastischer 
Verringerung des Schaltungsaufwandes, ihsbesondere der Verririgerung der Anzahl der bendtigten Speicherzel- 
len. 

Die Abtastfrequenz muB immer in Einklang mit dem augenblicklichen Frequenzinhalt — Abtasttheorem -r- 25 
des Filters stehen. Wird durch ein (digitales) TiefpaBjfilter der Frequenzinhalt gesenkt, so kann dementsprechend . 
auch die Abtastfrequenz hinter dem Filter reduziert werderi. Nichtrekursive Filter eigrien sich besonders gut fur 
diese Aufgabe,.weil sie in sehr einfacher Weise gestatten, die Ausgangstaktfrequenz zu sen ken: Man berechnet 
die nicht bendtigten Zwischenergebnisse einfach nichL Daraus ergibt sich zudem ein erheblicher Vorteil bezug- 
lich Aufwand und Geschwindigkeil So sind nur die Resultate des vorangegangenen Teilfilters bei entsprechen- 30 
der reduzierter Abtastrate von dem nachfolgenden Filter wet terzuverarbeiten. . 

Hier wird nun auf besondere Weise von der Abtastratenreduktion mit den Faktoren Rj bzw. n bei der 
Hintereinanderschaltung (Kaskadierung) yon niehreren Detektionsfiltern Gebrauch gemachL Die in Fig, 5a 
gezeigte Kaskadierungsstruktur eines N s -stufigen Mittelungsfilters mit seinen einzelnen Filterstufen der Lange 
Mi und deren diskreten StoBfunktionen M^k) ermoglicht eine besonders vorteilhafte Abtastratenstaffeluhg, ' 35 
wennfurdie jeweiligen Abtastraten-Reduktionsfaktorengesetzt wird: 


40 


45 


DJu die Abtastfrequenz f r am Ausgang des jeweiligen Filterabschnittes wird genau um den Faktor Mi des 
betroff enen Filterabschnittes abgesenkt und damit die nachfolgende Filterstufe getaktet. Dies zeigt sich beson- 
ders deutlich anhand der z-transfbrmierten Obertragungsfunktion M(z) des Gesamtfilters, Fig. 5b. Der Verzdge- 
rungsfaktor n — ist hier die negative Potenz der jeweiligen z-Variabein — erhalt man fiber die entsprechende 
TaktperiodendauerTi l/f r des entsprechenden Filterabschnittes nach folgender Beziehung: 50 

M i #T i 1 i 
r i+1 = — — = — H T n . (26) 

T T n==1 55 

Ebenfalls gilt hier n « 1. AuBerdem ist damit das Abtasttheorem immer genau erf Gilt, denn die reduzierten 
Filterbandbreiten, GL 10, der einzelnen Filterstufen sind immer etwas kleiner als die Halfte der Stufenabtastfre- 
quenzf r . 

Per besondere Vorteil dieses Verfahrens zeigt sich vor allem bei den stochastischen Filterparametern. So eo 
ergibt sich fur die (deterministische) Rechenwerksauflosung nur ein Wert von * 
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-N mit 


L = H M{ fur ±=1,2,3, ... ,N_ (27) 


wobei Lder Mittelungsproduktfaktordarstellt und mit derGesamtfilteriange identisch ist. Jedoch ist schon beim 
Effektivwert des Frequenzfehlers 


2 (i-N) jg> 
712 1 


= — — • — — mit 


(28) 


S = 2 S Mi 
i=l 


unter Verwendung des Mittelungssummenfaktors S eine wesentliche Verbesserung gegenuber den herkdmmli- 
25 chen Methoden,GL 12 usw^ gegeben. Gleichef maBen gilt dies auch fur den Signalstorabstand. . 



(29) 

S 

sowie fur die effektive Bitzahl 

D f —10 log Ii 

N f '» 1+ — (30) 

6,02 

wobei hier. fur Df der dB-Wert aus GI. 29 und fur L die skalare GroBe von GL 27 einzusetzen ist Daraus bestimmt 
sich dann die effektive Frequenzauflosung zii: 

q f = - — - • f a • 2 (3i) 


Folglich kann fur jedeh Einzelfall von kaskadierten Detektionsfiltern mit Abtastratenreduktion der Auflo- 
50 sungsgewinn ■ ' 

Aq f =2< N f- N ) 

gegenuber den Gmndverf ahren gleicher Filterlange, GL" 1 1 usw^ damit angegeben werden. 

55 Im folgenden wird nun eine Schaltungsanordnung eines derartigen Detektionsfilter anhand der Fig. 6 naher 
erlautert. Dabei handelt es sich urn ein vierstufiges Detektionsfilter fur 4mal 4bit-Phasenworte Oj(k) bis 04(k). 
Diese wurden aus dem analogen Eingangssignal s(t) beispieisweise in der beschriebenen Weise zunachst in 
einem 8-bit- Johnson-Code erzeugt, der anschlieBend in einen 4-bit-Binar-Code umgewandelt wurde. Insbeson- 
dere fur hohe und/oder sehr hohe Taktfrequenzen, z. B, groBer 100 MHz, ist es vorteilhafu die in serieller Form 

60 vorliegenden 4-bit-Binar-Phasenworte zunachst einem nicht dargestelltem 1 :4 Demultiplexer zuzuleiten, an 
dessen Ausgang dann in paralleler Form jeweils vier 4-bit-Binar-Phasenworte anliegen. Damit wird auf den 
Verbindungsleitungen von dem Ausgang des Demultiplexers zu dem Eingang (4>i(k) bis <I>4(k)) der Anordnung 
gemaB Fig. 6 die Obertragungsfrequenz vorteilhafterweise urn den Faktor 4 gesenkt, so daB schaltungstechni- 
sche Prdbleme, z. B. Obersprechen, vermieden werden. An dem Eingang des Filters liegen daher jeweils vier 

65 4-bit-Binar-Eingangsphasenworte <D(k) gleichzeitig an und werden dementsprechend auch gleichzeitig weiter- 
verarbeitet. Ohne den Verarbeitungsalgorithmus hierbei fur das Filter zu verandern, wird die Taktfrequenz 
bereits ab der ersten Stufe — Subtraktionsstufe — und fur alle nachfolgenden Stufen auf ein Viertel der 
ursprdnglichen Abtastfrequenz fA reduziert: Damit ist fur die fi it eri nterne Taktperiodend a uer T = 4/f A zu 


8 
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setzeri. Aus Grunden der Obersichtlichkeit wird fur die nachfolgenden Formeln jedoch vori der Originalabtast' 
frequenz fA bzw. vori der Origirialperiodendauer Ta » 1/Fa ausgegangen. In dem Beispiel ist eine Filterstruktur, 
mit N = 4 Filterstufen gewahh. Die dort mit T. 2T, 8T, 16T bezeichneten Bauelemente sind Verzogerungsglieder; ; 
weiche die anliegenden,Signale um die angegebene Anzahl von Taktperioden T=4-Ta verzogenvDas Bauele- 
ment T yerzogert daher mit T = 4-Ta, 2T verzogert mit 2T=2-4-Ta und so weiter. Hierbei stehen die Teilfilter- 
. langen: Mi « 8, M2 « 4, M3 = 2 und M4 = 2 durch ganzzahlige Potenzen der Zahl 2 im Verhaltnis. 
Dementsprechend gelten dann nach GL 25 fur die einzelnen Filterstufen die Abtastraten-Reduktionsfaktoren: R\ . 
= 1»R2 =* 8, R3 = 4 und R4 = 2 bzw. nach GL 26 fur die Abtastraten-Quotienten ( Verzogerungsfaktoren) n = 1, 
T2 = 8, n = 32, r4 = 64. Damit kann dahn die prinzipielle Rechenvorschrift zur Ermittelung des Frequenzwortes 
I - > ■ 

12.24 7 

' F(3c)= — • 2 2 2 2 L<f> (k— m^-m-j •ii^-m^) mit 

27T-128 ^4=1 m 3 =l m 2 =l m i =0 

A0(k) = [0(k-l)-0(k) ] mod 2tt (32) 

ohne jeglichen Zeitbezug angegeben werden. Dies reicht jedoch vollstandig aus, um die wichtigsten deterministic 20 - 
schen als auch stochastischen Detektbrparameter bestimmen zu kdnneri. So ergibt sich aus GL 27 eine Rechen- 
werksauflosung mit q 0 = 0,488 • 10~ 3 • fA- Hierbei besitzt der effektive Frequenzfehler nach GL28 den Wert 07 
= 1,127 • It)" 3 und die GL 29 liefert fur dieses Detektionsfilter einen Signalstorabstand von Df = 48.2 dB. Dies 
hat nach GL 30 eine effektive Bitzahl Nf = 5J5 bit zur Folge, mit der sich nach GL 31 die effektive Frequenzauf Id- 
sung zu qr » 0,173 • 10"? - fA bestimmen laBt Damit besitzt dieses Filter eine um das 2 f 8fache bessere* effektive 25 
Frequenzauflosung als es sich aus der Rechenwerksauflosungq 0 direkt ergibt. 

. , Dagegen lassen sich die zeit- und frequenzabhangigen Eigenschaf ten des DetektionsfHters am besten anhand 
seiner z.-transformierten Obertragungsfunktionen 
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15 


M(z) « M 1 (z).M 2 (2).M 3 (z).M 4 (z) mit 
H x (z) = — (2°+z" 1 +z r2 +z" 3 +2" 4 +z" 5 +z" 6 +z" 7 ) 7 
M 2 (z) = -— z^ 7 (z°+z" 8 +z" 16 fz- 24 ) , (33) 
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1 


M 3 (z) - — - z r31 (z°+z" 32 ) und 


M 4 (z) = z^ 3 (z 0 4-z" 64 ) 


50 


fur das Gesamtfilter M(z) sowie fur die Filterstufen Mi(z) erklaren. Fur die Variable z wird hierbei die ailgemein . 
bekannte Grundbeziehung z = exp Q2kt/T) der Z-Transformation verwendet. Zur Herleitung der diskreten und 
auf fA « 1 Hz normierten Amplitudengange G*(K) — |Mi(j27ik- 1)| aus GL 33 fur das Filter und deren Teilfiltern 
sind die Formeln der Euler'schen Relation fur z n ± z~ n zweckmaBig. Fur eine gerade Filterkoeffizientenanzahl 55 
Mi ergibt sich somit ein Amplitudengang: " . 


1 

Gi (k)=2- 7. cos [t-k-riCn-rO, 5) ] mit 1^/2 (34) 
und fur eine ungerade Filterkoeffizientenanzahl' Mi: 
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V G f (k)=l+2. 2 X cosCt-k-ri-n] mit ^-(M-l)/^ . (35) 
n=l 


Unter Verwenduhg der Gl. 34 und 35 werden nun in Fig. 7 fur dieses 4-stufige Detektionsfilter die normierten 
Amplitudengange fur a) mit Abtastraten-Reduktionen und b) ohne Abtastraten Reduktion beziiglich der.redu- 
zierten und normierten Abtastbandbreite f r - 1/128 des FrequenzmeBwertes F(k) aufgezeigt Die in den 
Kopfzeilen aufgeiisteten Parameter beziehen sich hierbei auf das Filter mit Abtastraten-Reduktion. 

Bei Vergleich der beiden Kurven in Fig. 7 stellt sich vor allem ein erheblicher Unterschied in den 3dB-Band- 
breiten heraus. Die normierte Bandbreite des Filters mit Abtastratenreduktion betragt hierbei Fn = 0,000211 
und 1st sbmit zur Bandbreite. GL 10, des aquivalenten Filters ohne Abtastraten-Reduktion um den Wert 18,5mal 
kleiner. Dieser groBen Verbesserung in der Frequenzauflosungs-Bandbrehe steht dem nur eine relativ geringe 
Zunahme der MeBzeit von 127 -T A auf 164,5-Ta entgegen. Jedoch ist es in der technischen Realisferung leicht 
moglich, durch ein optimales Timing" des Rechenablaufes im Filter diese MeBzeitzunahme von ca. 36 Abtast- 
taktperioden T A erheblich zu verringern. Damit erbringt die Kaskadenlosung mit entsprechend reduzierten 
Abtastratenfaktoren, Gl. 25 und 26, fur die einzelnen Filterstufen eine wesentliche Steigerung der (statistischen) 
MeBgenauigkeh ohne jedoch nennenswerte Verschlechterungen bei der MeBzeit zu erhalten. AuBerdem ergibt 
dieses Verfahren bei der tephnischen Ausfuhrungden geringsten Aufwand an Spetcherzellen und Addierstufen. 

Die Erfindung ist nicht auf die beschriebenen Beispiele beschraiikt, sondern sinngemaB 'auf weitere anwend- 
bar. Beispielsweise ermoglichen digitale GaAs-Schaltkreise, monolithisch integrierte Mikrowellen-Bausteine 
(MMIC) und/oder ECL-Schaltkreise im Picosekundenbereich (Gijgabit-Logik) in steigendem MaB eine digitale 
Signalverarbeituhg von hohen und sehr hohen Frequenzen (>0,5 GHz). Damit ergeben sich fur das hier 
beschriebene Momentanfrequenz-Detektionsverfahren in Verbindung mit der Phasen/Digital-Umsetzung im- 
mer breitere und heue Anwendungsmoglichkeiten auf alien Gebieten der Kommunikatiohs : und Radartechnik. 
Hinzu kommt noch die erhebliche Verringerurig des technischen Aufwandes bei der Realisierung dieses Detek- 
/rionsfilter-Prinzips. So entfallen unter anderem jegliche Koeffizienten-Multiplizierer, da nur Filterkoeffizienteri 
mit —1,0 und + 1 auftreten konnen. Somit ermoglicht es eine einfache und extrem schnelle Verarbeitung yon 
quantisierterTPhasenwinkelwerten mit reellen Algorithmen, die besonders gunstig in Hochgeschwindtgkeits- 
ECL- oder GaAs-Gate- Arrays implementiert werden konnen. 

Patentanspriiche 

1. Verfahren zur Momentanfrequenz-Detektion, bei dem aus einem analogem Eingangssignal mittels eines 
Taktsignales digitale Abtastwerte gebildet werden und aus den Abtastwerten ein digitales Frequenzwbrt, 
das einem Frequenzanteil des Eingangssighales entspricht, ermittel t wird, dadurch gekennzeichne 

— daB durch eine quantisierte Phasenabtastung mittels eines Taktsignales zunachst zu jedem Takt ein 
Phasenwort gebildet wird, 

— daB aus jedem Phasenwort mittels eines 2errKompIement-Kodierers zwei als 2er-Komplement 
. kodierte Phasenworte gebildet werden, 

— daB aus den im 2er-Komplement kodierten Phasenworten durch eine Verzogerungsschaltung 
verzogerte und unverzogerte kodierte Phasenworte erzeugt werden, 

— daB durch eine digitale Addition der verzogerten und der unverzogerten kodierten Phasenworte 
zunachst vprlaufige Frequenzworte erzeugt werden und 

— daB die vorlaufigen Frequenzworte mittels einer digitalen Mittelwe^rtsbildung in endgOltige Fre- 
quenzworte umgewandelt werden. 

2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daB die Mittelwertsbildung ausschlieBlich mit Hilfe 
von Verzogerungs- und Additionsstufen erfolgt 

3. Verfahren nach Anspruch 1 oder Anspruch 2, dadurch gekennzeichnet, daB aus den Verzogerungs- urid 
Additionsstufen ein digitales Filter mit kanonischer Struktur gebildet wird. 

4. Verfahren nach Anspruch 1 oder Anspruch 2, dadurch gekennzeichnet, daB aus den Verzogerungs- und 
Additionsstufen ein digitales Filter mit nicht kanonischer Struktur gebildet wird 

5. Verfahren nach einem der vorhergehenden Anspruche, dadurch gekennzeichnet, daB aus den Verzoge- 
rungs- und Additionsstufen ein digitales Filter, das eine Blocksummenaddition en thai t, gebildet wird. 

6. Verfahren nach einem der vorhergehenden Anspruche, dadurch gekennzeichnet, daB die Mittelwertsbil- 
dung mittels einer Kaskadierung von mindestens zwei digitalen Filtern, die jeweils als TiefpaBfilter ausge- 
bildet werden, vorgenommen wird. 

7. Verfahren nach Anspruch 6, dadurch gekennzeichnet, daB bei der Kaskadierung nach jedem digitalem 
Filter eine Verringerung der zum Betreiben der Filter erforderlichen Taktrate erfolgt 

8. Verfahren nach einem der vorhergehenden Anspruche, dadurch gekennzeichnet, daB bei einer Addition 
zweier im 2er-Komplement vorliegender Phasenworte ein Obertrags-Bit (Carry-Bit) hinzugefugt wird. 

9. Verfahren nach einem der vorhergehenden Anspruche, dadurch gekennzeichnet, daB die Mittelwertsbil- 
dung mit in paralleler Form vorliegenden Phasenworten erfolgt 

10. Verfahren nach einem der vorhergehenden Anspruche, dadurch gekennzeichnet, daB bei der Auswer- 
tung eines als BandpaBsigna! vorliegendem analogem Eingangssignal zu den endgultigeh digitalen Fre- 
quenzworten ein digitaler Frequenz-Korrekturwert welcher die Frequenzlage des BandpaBsignales bein- 
haltet, addiert wird. 
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» . 11. Verfahren nach einem der vorhergehenden Ansprflche; ;zii Venyendiirig bei der Detektion und/dder 

Demodulation von frequenzmodulierten und/oderphasenmoduliertenanalogen Eingangssignalen; 
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